Anisotropic rare-earth spin ensemble strongly coupled to a superconducting resonator 
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Interfacing photonic and solid-state qubits within a hybrid quantum architecture offers a promis- 
ing route towards large scale distributed quantum computing. Ideal candidates for coherent qubit 
interconversion are optically active spins magnetically coupled to a superconducting resonator. We 
report on a cavity QED experiment with magnetically anisotropic Er^+:Y2Si05 crystals and demon- 
strate strong coupling of rare-earth spins to a lumped element resonator. In addition, the electron 
spin resonance and relaxation dynamics of the erbium spins are detected via direct microwave ab- 
sorption, without aid of a cavity. 



Quantum communication networks are considered to 
distribute entangled states over a large scale computing 
architecture [TJ |2]. The core elements of future quan- 
tum networks, i.e. quantum repeaters [3] as well as net- 
work nodes, will be realized by using qubits and quan- 
tum memories of diverse physical nature [H [5] . Today, 
elementary quantum networks linking two remote single 
atoms have been demonstrated [H [?]• Solid state sys- 
tems, such as superconducting (SC) quantum circuits [8], 
nanomechanical devices [9 and spin doped solids [10] po- 
tentially offer larger scalability and faster operation time 
compared to systems based on the single atom approach. 
However, such solid state devices operate at microwave 
and RF frequencies, which are less suitable for long-range 
quantum communication than optical channels due to 
losses in cables and the high noise temperature of anten- 
nas. Therefore, to establish a fiber-optical link between 
them, one has to use a quantum media converter, i.e. a 
device which coherently interfaces matter and photonic 
qubits [m© . 

One of the promising ways towards implementation of 
such a converter relies on using optically active spin en- 
sembles in a hybrid quantum architecture p!5HT7] . At the 
present time, the research activity is primarily focused 
on various cavity QED experiments with spin ensembles 
of NV-centers in diamond p!8H22j . Recently, collective 
strong coupling has been demonstrated in conventional 
electron spin resonance (ESR) experiments with organic 
molecules in a 3D cavity [23, 24 , and four- wave mixing 
from Fe^+ ions in a 3D sapphire loaded cavity has been 
shown [25]. Surprisingly, less progress has been achieved 
with paramagnetic laser materials such as ruby [19 and 
rare-earth ions doped crystals [T7l|26]. The way to the 
strong coupling regime in these initial experiments was 
impeded by the large inhomogeneous broadening of spin 
ensembles [27] . 

Unlike electronic spins of 3(i transition series metal 
ions, Kramers rare-earth (RE) ions reveal a rather strong 
magnetic anisotropy due to the distortion of their 4/ 
electronic orbitals by a crystal field, see Chapter 5 of 
Ref. |28^. Such a distortion induces a strong dependence 



of the g-factors and the Rabi-frequency on the orientation 
of the RE-ion doped crystal with respect to the polar- 
izations of the permanent (DC) and the oscillating (AC) 
magnetic fields [29] . The magnetic moments of Er^+ ions 
in certain crystals with axial symmetry can vary from 
zero to 7.5/iB, where /i^ is the Bohr magneton [30} [3T]. 
Using such magnetically anisotropic materials in hybrid 
quantum systems leads to new challenges. In this let- 
ter, we present cavity QED experiments on anisotropic 
Er^+:Y2Si05 (Er:YSO) crystals and demonstrate strong 
collective coupling of the RE-spin ensemble to a super- 
conducting lumped element (LE) resonator. 

The experimental setup is shown in Fig. [ija). Two 
Er:YSO crystals are glued on a superconducting niobium 
chip containing 9 lumped element resonators coupled to a 
50 1^ transmission line [32] . The resonance frequencies of 
the circuit occupy the microwave band between 4.5 and 
5.2 GHz and exhibit loaded and intrinsic quality factors 
of Qc — 10^ and Qi 10^, respectively. The quality 
factors of the LE resonators do not degrade significantly 
by applying an in-plane DC magnetic field up to 340 mT. 
Figure [ijb) presents a simulation of the AC magnetic 
field and the picture of a single LE resonator. Due to the 
narrow width of 10 /im of the meandered inductance, the 
AC microwave field Bi cos ut is mainly concentrated in 
direct vicinity of the chip. The mode volume of such a 
resonator is Vm ^ 10~^ cm^ comprising Ng ~ 10^^ spins 
for magnetic coupling. 

We study two Y2Si05 crystals (noted below as 
EriYSOl and Er:YS02) doped with 200 ppm Et^+ ions 
(Scientific Materials, Inc.). The orientation of their 
optical extinction axes 6, Di and D2 is sketched in 
Fig. [TJc) [31 . The alignment of the crystals in the 
DC magnetic field Bq is specified by the angles 6 and 
(j). The comparison of measured ESR spectra to the 
"Easyspin" simulation [33] yields the following pair of 
angles 61 = 44°, = 111° and O2 = 48°, (/)2 = 85° for 
the Er:YS01 and Er:YS02 crystals, respectively. 

Figure [TJd) qualitatively outlines the main idea of the 
experiment: the EriYSO crystal is known to have the 
strongest magnetic anisotropy due to its low axial sym- 
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metry C2h^ and the principal values of the g-tensor for the 
crystallographic site 1 are ^ 0, g^ = 1.5 and g^ = 14.8, 
see [30]. The angular dependence of the DC and AC g- 
factors (g and g^^) of the crystal rotated around x-axis of 
the g-tensor are given by (s.f. Chapter 3 of Ref. ^28^ ) 
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The coupling strength of a single spin to a microwave 
photon with an AC field of Bi cos ujt is determined by 
'^1 = l~^Bgi{^)\Bi\/2fi. Therefore, to attain the largest 
coupling strength, the AC field has to be polarized along 
the g^ direction and the DC field should be aligned per- 
pendicular to the AC field. In the following we show, 
that in order to reach the strong coupling regime one has 
to sacrifice the high spin tuning rate of EriYSO. 

The experimental system is placed inside a BlueFors 
LD-250 dilution fridge and the on-chip ESR spectroscopy 
is performed in a temperature range of 20-500 mK. Fig- 
ure |2] presents the transmission spectrum of the circuit 
|5'2i(cj)| as a function of the applied magnetic field. The 
excitation power at the input of the circuit is approxi- 
mately 50 aW corresponding to a few microwave photon 
excitation level inside each resonator. Every horizontal 
line in the ESR spectrum corresponds to one of the 9 
LE resonators on the chip. Below 200 mT, the reso- 
nance lines are interrupted by a regular pattern of dis- 
persive cavity shifts due to the coupling to the magnetic 
dipole transitions of Er^+ between the electronic states 
rris ±1/2. 

In general, the ESR spectrum of EriYSO consists of 
two pairs of strong lines associated with the occupation of 
ions of two distinct crystallographic sites and each 
in two magnetically inequivalent positions [30l [311 [34] . 
The dashed and dotted lines in Fig. [2] correspond to the 
spin tuning lines of the samples EriYSOl and Er:YS02, 
respectively. The term S'l^, for instance denotes the 
transition of the erbium ions in site 1 and magnetic in- 
equivalent position (magnetic class) a. Correspondingly, 
the term S2i) belongs to site 2 and magnetic class h. The 
effect of the magnetic anisotropy is clearly seen in the 
spectrum in Fig. [2] High field transitions with smaller 
g-factor have larger coupling strength. 

Inset (a) in Fig. [2] shows the positions of the resonators 
beneath the crystals which were partially identified by 
measuring the local photoresponse of the circuit in a 
laser-scanning- microscope (LSM) setup [35]. A typical 
LSM refiectivity map taken at the excitation frequency 
of resonator #5 is presented in inset (b). These reflectiv- 
ity maps have been taken during a separate run at 4 K, 
and only an incomplete identification of the resonators 
was possible. 

In the field region between 210 and 265 mT, a faint 
narrow absorption line traverses the spectrum, which cor- 
responds to the S2a transition of the Er:YS02 crystal. 



The interception of this line with resonances #6-8 results 
in a complex hybridization pattern, due to the mutual in- 
ductive coupling between those resonators. 

Figure [3] presents the detailed study of the avoided 
level crossing between the erbium transition S2a and res- 
onator #5 at 250.6 mT. The raw transmission spectrum 
\S2i{^)\^ of the circuit measured at the avoided level 
crossing at 253.1 mT is shown in Fig. [sjb). For com- 
parison, we show the transmission spectrum taken away 
from the anticrossing at 204.3 mT. A typical experimen- 
tal power spectrum \S2i{oj)\^ measured at the center of 
the avoided level crossing consists of the normal mode 
splitting itself, 2 additional uncoupled resonances and the 
absorption dip. The spectrum \S2i{oo)\^ can be described 
by the following expression 
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where B{uo) is a second order polynomial, accounting for 
the baseline, uoi are the frequencies of the resonances and 
7^ are their linewidths. The coefficients an and a2i de- 
scribe the absorptive and the dispersive parts of the res- 
onances. The fit of the experimental data recorded at 
251.3 mT to Eq.|3]is presented in Fig.jsjb). 

The dips marked with |+) and |— ) correspond to the 
normal mode splitting. The #3 and #4 denote the posi- 
tions of the resonators which are not coupled to Er:YS02. 
In order to extract the spectrum of the normal mode 
splitting we normalize the measured curve \S2i{oj)\^ to 
the baseline B{uo) and subtract the non- interacting res- 
onators #3 and #4 as well as the transmission line ab- 
sorption dip. Figure [sjc) shows the corrected spectra 
measured at field values of 251.3 mT (gray line) and 
at 250.6 mT (dark gray line). Two well-separated dips 
are confirming strong coupling of the spin ensemble S2a 
to the cavity. To describe the normal mode splitting 
we use the model of coupled quantum harmonic oscil- 
lators [ini [iSl [36]. The corrected spectra are fit to 
the avoided level crossing curves according to the equa- 
tions given in Ref. [19 and presented in Fig. [sj^c) by 
dashed lines. The external kc/2TT = 4.7 MHz and in- 
ternal ki /27r = 0.7 MHz damping rates of resonator #5 
were determined away from the avoided level crossing at 
a field of 204.3 mT and set constant for the fit. We ob- 
tain a splitting size of 2vl2i\ = 68±1 MHz and a spin 
linewidth (FWHM) of F^/tt = 24±0.5 MHz, which gives 
the cooperativity parameter C = 2v'^ / KcV^ ~ 36. 

The inset in Fig. [sj^a) shows the infiuence of the 
i67g^3+ Jiyperfine transitions on the damping rate of 
resonator #5. Particularly at frequencies below 6 GHz 
the hyperfine spectrum is rather complex and consists 
of about 20 lines in the field range of our interest (205- 
280 mT). We could not identify all magnetic hyperfine 
transitions of EriYSO with sufficient confidence. Above 
6 GHz the hyperfine spectrum appears as a regular pat- 
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tern of lines and the classification of the spin transitions 
is easier [171 130] • Nonetheless, the transitions between 
states with equal nuclear spin projection Am/ = (HF 
in the inset) possess larger coupling strengths than 
quadrupole transitions with Am/ = ±1 (Q in the inset). 
The fit of the damping rate of the cavity /ci/27r, weakly 
coupled to HFl transition, yields v^/i/27r = 4 MHz cou- 
pling and a linewidth of r^^i/27r = 7.6 MHz (see Ref. [17] 
for the fitting procedure). For the quadrupole transition 
Q the coupling strength of spins and their linewidth are 
Vq/27i =1.8 MHz and r^/27r = 7.1 MHz respectively. 
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TABLE I: g-factors, coupling strengths v and linewidths r2 in 
MHz of 4 spin transitions of the Er:YS01 coupled to resonator 
#9 (left), and Er:YS02 coupled to resonator #5 (right). 



Low-field electronic spin transitions of the EriYSOl 
and the Er:YS02 crystals do not reveal the strong cou- 
pling; their coupling strengths and linewidths are found 
from dispersive shift of the cavity. The results are sum- 
marized in the Table |lj The coupling strength is in- 
versely proportional to the g-factor. The spin linewidth 
for the S2a transition for Er:YS02 drops to 12 MHz 
after an initial increase. The study of Er:YS02 crys- 
tal on the X-band ESR spectrometer Bruker Elexsys 
580 at a temperature of 6.3 K suggests a minimal 
linewidth r2/27r o:^ 14 MHz for low field transitions and 
r2/27r 25 MHz for the high field ones. The measured 
behavior of the spin linewidth in the on-chip ESR exper- 
iment can be explained by the spatial inhomogeneity of 
the magnetic field in the vicinity of the SC chip due to 
the Meissner effect. The angular variations of the mag- 
netic field broadens mainly those transitions which are 
close to the sharp vertex of the g-factor ellipse, see also 
Fig-Bd). 

The measured narrow linewidth of the low-field erbium 
transition S2a is consistent with the measurements of the 
on-chip microwave absorption. In the absence of any res- 
onator, the microwave signal is extinguished just due to 
the proximity of the spins to the transmission line. Such 
an effect has recently been demonstrated with a highly 
doped ruby crystal [19]. In contrast, the spin concentra- 
tion in our experiment is UEr — 7x 10^^ cm~^, which is at 
least by a factor of 10 lower. A clear signal results from 
microwave absorption by Ng ~ 10^^ spins with large AC 
g-factor of nearly 15. 

We studied the absorption line at a magnetic field 
of 273.2 mT, which corresponds to an ESR frequency 
of 5.331 GHz (440 MHz away from resonator #5). In 
order to distill the absorption profile from the trans- 



mitted signal |S'2iP, we normalized the latter to the 
far detuned (140 MHz) transmission taken at 266 mT. 
The resulting absorption spectrum recorded at an ex- 
citation power of 0.1 fW is shown in Fig. [4]^ a). The 
experimental data (gray dots) show 12% absorption 
of the transmitted microwave signal and fit well to 
a Lorentzian (solid line) with an ESR inhomogeneous 
linewidth r^/27r = 13.8±0.4 MHz. The linewidth is 
slightly increased compared to the avoided level cross- 
ing due to the larger geometrical size of the transmission 
line, which makes it more sensitive to the magnetic field 
inhomogeneity. Therefore, a linewidth narrowing due to 
the cavity protection effect can be excluded here [37 . 

The relaxation time of the spin ensemble at mK tem- 
peratures has been recently measured for NV-centers by 
using the dispersive cavity shift [20 . In our experiment, 
the large absorption signal allows for observing the spin 
relaxation dynamics directly. For that purpose the vec- 
tor network analyzer was tuned to the absorption dip 
at uja/27r = 5.331 GHz. The probing power was set to 
60 fW and the chip was irradiated for 1 sec with an in- 
tense 0.1 nW pulse at a central frequency of 5.330 MHz. 
The reappearance of the absorption signal is presented in 
Fig. ^h). The magnitude of the absorption dip is pro- 
portional to the average population difference of the spin 
ensemble, and therefore, to its total magnetization 

1^21 1'(0) - 1^21 CX iVi - iV2 cx M,{t). (4) 

After the pulse, the population restores its equilib- 
rium level N2/N1 = exp{-hLjJa/kBT) ^ 10"^, where 
T = 20 mK is the temperature of the experiment. The 
solid line in Fig. ^h) shows the fit of the measured 
spin relaxation signal to the exponential decay with 
Ti = 4.3±0.2 sec. 

To conclude, we have presented cavity QED experi- 
ments with magnetically anisotropic Er^+:Y2Si05 crys- 
tals. The coupling strength and linewidth strongly de- 
pend on the g-factors of the erbium spin transitions. The 
narrow inhomogeneous linewidth of 12 MHz for the high 
field transition with g = 1.4 allowed us to attain strong 
coupling between the SC resonator and the RE-spin en- 
semble. The ESR as well as the spin relaxation dynamics 
were detected by direct measurement of the microwave 
absorption of the spins coupled exclusively to the trans- 
mission line. The presented experiment demonstrates the 
promising potential of rare-earth ion doped crystals for 
application in hybrid quantum architectures. 
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Rotation angle, cp 



FIG. 1: (Color online) (a) Picture of the experiment: two Er:YSO crystals are placed on the SC-chip with 9 LE resonators. 
The DC magnetic field is applied along the chip's surface, (b) The simulated AC magnetic field in the plane perpendicular to 
the DC field in the vicinity of the LE inductor (meander structure). (Inset) Picture of a LE resonator, (c) Orientation of the 
crystal's axis, (d) Dependence of g and g^^ on the rotation angle of the crystal around the x-axis of the g-tensor. The maximal 
coupling is reached when the AC field is aligned along the largest component of the g-tensor. 





FIG. 2: (Color online) ESR spectrum of the Er:YSO crystals coupled to the LE chip. The dashed lines correspond to 4 
electronic spin transitions of the Er:YS01 crystal and the dotted lines to transitions of the Er:YS02 crystal. The resonator 
labels are given on the right. (Inset a) Positions of the LE resonators beneath the crystals. (Inset b) LSM reflectivity map 
measured at the excitation frequency of resonator #5. 
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FIG. 3: (Color online) (a) Transmission spectrum of resonator #5 strongly coupled to the electronic spin S2a transition. 
(Inset) The internal damping rate of the resonator ki /27r clearly resolves the hyperfine spectrum of ^^^Er"^^. The "HF's" 
denote magnetic hyperfine transitions and "Q" marks a quadrupole transition, (b) The gray dashed line is the transmitted 
power P21 = |*S'2i(ct;)|^ measured away from the avoided level crossing at 204.3 mT. The dark gray line displays the power 
spectrum at the anticrossing at 251.3 mT, which is fit using Eq. [3] (blue line), (c) The corrected transmitted power spectra 
taken at 250.6 mT (dark gray line) and at 251.3 mT (gray line) clearly show a normal mode splitting. The dashed lines show 
the fit to this splitting. 
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FIG. 4: (Color online) (a) The gray curve presents the absorption profile of the erbium spins coupled to the transmission line 
at 273.2 mT. The solid line shows a Lorentzian fit to the data, (b) Reappearance of the absorption signal (gray curve) at 
5.331 GHz after irradiation of the chip with an intense microwave pulse. The solid line is a fit of the data to an exponential 
decay. 



